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Experimentally detected ultrafast spin-avalanches spreading in crystals of molecular 
(nano)magnets (Decelle et al., Phys. Rev. Lett. 102 , 027203 (2009)), Ref. [T], have been re¬ 
cently explained in terms of magnetic detonation (Modestov et al., Phys. Rev. Lett. 107 , 207208 
(2011)), Ref. [2|. Here magnetic detonation structure is investigated by taking into account trans¬ 
port processes of the crystals such as thermal conduction and volume viscosity. In contrast to the 
previously suggested model, the transport processes result in smooth profiles of the most impor¬ 
tant thermodynamical crystal parameters - such as temperature, density and pressure - all over 
the magnetic detonation front including the leading shock, which is one of the key regions of mag¬ 
netic detonation. In the case of zero volume viscosity, thermal conduction leads to an isothermal 
discontinuity instead of the shock, for which temperature is continuous while density and pressure 
experience jump. 

PACS numbers: 


I. INTRODUCTION 

Presently there is much interest in molecular (nano) 
magnets with unique superparamagnetic properties, 
which may be used for quantum computing and memory 
storage! 3 7 A remarkable feature of nanomagnets is that, 
in contrast to classical magnets, these macromolecules 
with large effective molecular spin (e.g., S = 10 for Mni 2 - 
acetate) can keep their spin-orientation upon the reversal 
of the external magnetic fieldBecause of the strong 
molecular anisotropy, spin of a nanomagnet is directed 
preferentially along the so-called easy axis, and it leads 
to a considerable energy barrier between the spin-up and 
spin-down states. At low temperatures, in a magnetic 
field directed along the easy axis, the states with spin 
along the field and against the field become stable and 
metastable, respectively. The energy difference between 
the two states is determined by the Zeeman energy, Q , 
as illustrated in Fig. |T| with the energy barrier desig¬ 
nated by E a . The barrier hinders spontaneous quan¬ 
tum tunneling from the metastable to stable state at low 
temperatures,®" 3 ^ so that fast spin-flipping requires help 
from outside. 

For nanomagnets composing a crystal, relatively fast 
spin-flipping of one particular molecule may be induced 
by energy supplied by its’ neighbors. When all or most of 
the molecules of a crystal are initially in the metastable 
state, then local heating by an external source may trig¬ 
ger local spin-flipping, with Zeeman energy released in 
the heated r egion and transported to the next layer of 
the crystal EHnUEU The heat facilitates spin-flipping in 
the next layer, and so on, so that the process spreads in 
a crystal as a thin self-supporting magnetization front - 
spin-avalanche - well-localized spatially. Usually, energy 
in spin-avalanches is transported from one crystal layer 
to another by means of thermal conduction, and hence a 
spin-avalanche propagates at moderate speed, ~ 1 — 10 
m/sPHZ Due to striking similarity of such avalanches to 
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FIG. 1: (Color online) (a) Schematic representation of the 
double-well structure of a nanomagnet, (b) Energy levels for 
Mni 2 -acetate in the external magnetic field B z = IT. Axis z 
is parallel to the easy axis of the crystal. The energy barrier 
(activation energy) E a and the Zeeman energy Q are indi¬ 
cated. 


slow combustion flame, deflagration , avalanches of this 
type are typically called ’’magnetic deflagration”. 

In contrast to the slow magnetic deflagration s tudi ed 
in the absolute majority of works on the subject EHZmHH] 
recent experiments of Ref. [TJ detected ultrafast spin- 
avalanches propagating at speed comparable to the sound 
speed in the crystals, ~ 2000 m/s. The theory presented 
in Ref. [2] has explained the ultrafast spin-avalanches 
in terms of ’’magnetic detonation ” and investigated the 
key properties of the process. In particular, it has been 
demonstrated that magnetic detonation belongs to the 
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type of weak detonations and propagates with speed only 
slightly exceeding the sound speed « 2000 m/s. Paper [2] 
has also studied structure of magnetic detonation within 
the traditional combustion model of a detonation front 
consisting of an infinitely thin leading shock and a zone 
of energy release of finite thickness.® Such a model has 
been originally developed for gaseous detonations, which 
are quite strong and propagate with speed exceeding the 
sound speed in gases by an order of magnitude. However, 
it is not rigorous applying such a model to magnetic det¬ 
onation, which is extremely weak by combustion scales 
(although quite strong when compared to magnetic defla¬ 
gration). In contrast to strong shocks in gases, which are 
infinitely thin from the hydrodynamic point of view, a 
weak shock exhibits a continuous structure controlled by 
transportprocesses such as thermal conduction and/or 
viscosity.®^ The same property should be naturally ex¬ 
pected for magnetic detonation. 

The purpose of the present work is to provide accurate 
description of magnetic detonation structure in crystals 
of nanomagnets by taking into account thermal conduc¬ 
tion and volume viscosity. Here we show that, in con¬ 
trast to the previously suggested model of Ref. [2], the 
transport processes result in smooth profiles of the most 
important thermodynamical crystal parameters - such as 
temperature, density and pressure - all over the magnetic 
detonation front including the leading weak shock, which 
is one of the key regions of magnetic detonation. In the 
case of zero volume viscosity, however, thermal conduc¬ 
tion leads to an isothermal discontinuity instead of the 
shock, for which temperature is continuous while density 
and pressure experience jump. 

The present paper is organized as follows. In Sec. II we 
overview basic features of gaseous combustion detonation 
needed for proper understanding of magnetic detonation 
in crystals of nanomagnets. In Sec. Ill we present basic 
equations describing magnetic detonation, develop the 
analytical theory for the most important detonation pa¬ 
rameters and study magnetic detonation structure con¬ 
trolled by thermal conduction only assuming zero volume 
viscosity - such a structure involves isothermal disconti¬ 
nuity instead of a shock employed in Ref. [2J. In Sec. IV 
we demonstrate dramatic modifications of the magnetic 
detonation structure due to volume viscosity. 


II. BASIC FEATURES OF GASEOUS 
COMBUSTION DETONATION 


In this section we remind briefly the most important 
features and the methods of investigation of traditional 
gaseous combustion detonations, in order to highlight 
the similarity and difference between the combustion and 
magnetic detonations. By definition, detonation is a fast 
supersonic combustion regime, for which preheating of 
the cold fuel mixture happens due to the leading shock, 
see Fig. [2] In turn, shock propagation without decay in 
combustion detonation is supported by energy release in 
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FIG. 2: Schematic of strong gaseous detonation with char¬ 
acteristic profiles of temperature and pressure. 


chemical reactions in the active reaction zone. Typically 
for combustion, the activation energy of the reactions is 
quite high, so that reactions develop relatively slow just 
after the leading shock in the so-called induction zone 
(still, much faster than in the fuel mixture). It requires 
certain (induction) time for self-acceleration of chemical 
reactions in a gas parcel after passing the shock, after 
which the chemical reactions become fast and convert 
the fuel mixture into the burning products with release 
of a large amount of energy and strong expansion of the 
burning gas. 

Here we are interested in a planar stationary one- 
dimensional detonation propagating with constant super¬ 
sonic speed D in a uniform gaseous mixture; we take the 
detonation front propagating along the z-axis in the neg¬ 
ative direction. By adopting the reference frame of the 
stationary detonation front, we obtain the fuel mixture 
moving with velocity uq = D towards the leading shock 
of zero thickness; the initial density, pressure and temper¬ 
ature are designated by po , Pq and Tip The fuel mixture 
is compressed in the shock (label ”s”) with strong in¬ 
crease of density, pressure and temperature to p s > po, 
P s > Tq, T s > Ty, and with drop of the gas velocity to 
a new subsonic value u s < D. Temperature increase ini¬ 
tiates combustion reactions with release of chemical en¬ 
ergy, hence leading to even stronger temperature increase 
until it reaches the final maximal value in the detonation 
products (label ”d”), Td > T s . At the same time, pres¬ 
sure and density decrease in the reaction region, pd < p s , 
Pd < P 8 , as illustrated schematically on Fig. [2j 

Modifications of the gas parameters in a detonation 
front are described by the hydrodynamic laws of mass, 
momentum and energy conservation expressed in terms 
of fluxeJ23 


pu = consti = p 0 D, 


(1) 


P + pu 2 = const 2 = Pq + PoD 2 , 


( 2 ) 
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pu [ h + — I = constz = pqD h 0 + 


D 2 


where h is the gas enthalpy 

7 P 


h = 


( 7-1 )P 


+ aQ , 


(3) 


(4) 


7 is the adiabatic exponent, Q is the chemical energy 
stored per unit mass, a is fraction of the unburned fuel 
mixture, which changes from ao = a s = 1 in the fresh and 
shocked gas to a d = 0 in the detonation products. The 
hydrodynamic equations Eqs. 0 © are complemented 
by the equation of state, which is taken for gaseous det¬ 
onations in the form of the ideal gas law 

P=—RT, (5) 

m 

where R is the ideal gas constant and m is the molar 
mass of the gas (for simplicity we assume here that the 
molar mass does not change in the combustion process). 
Traditionally, the theory of shock waves and detonations 
employs volume per unit mass V = 1/p for analyzing the 

process P3 

Modifications of the gas parameters in a leading shock 
are described by Eqs. 0-0 with zero energy release 
in the reaction within the discontinuous shock, a s = 1. 
Then Eqs. 0-0 may be reduced to the so-called Hugo- 
niot equation, which specifies all possible finite states of 
the shocked gas P s , V s for a fixed initial state Po, Vb; the 
result is demonstrated by the blue solid curve in Fig. [3] 
for the initial state corresponding to the ideal gas at ini¬ 
tial volume per unit mass Vo = 6.33 m 3 /kg and pressure 
Po = 1.33 • 10 4 Pa. For any particular final state s at 
the shock, equations 0, (p]) relate the slope of the green 
straight line ( sO ) in Fig. plto the front propagation speed 
D asP 

** - fcv.- (6 » 

In order to obtain the final state of the detonation 
products we solve Eqs. 0-0 for the case of complete 
burning a d = 0 ; the result is plotted by the solid red curve 
in Fig. [3]for the scaled energy release ( 7 —l)Q/( 7 PoVo) = 
9 corresponding to the acetylene-air combustionP The 
red curve (for detonation) corresponds to higher pressure 
for the same volume than the blue curve (for a shock) due 
to the energy release in a detonation front. Still, unlike 
a shock, detonation speed is not a free parameter, but 
it is determined by particular boundary conditions. For 
example, the case of a freely propagating detonation cor¬ 
responds to the so-called Chapman-Jouguet (CJ) detona¬ 
tion regime, with the detonation products moving away 
from the front at local sound speed ! 21 The CJ detonation 
propagates with the smallest speed possible for detona¬ 
tions for a fixed energy release, and indicated by the green 
tangent line in Fig. [3] Thus, on the pressure-volume di¬ 
agram, the CJ detonation structure corresponds first to 



FIG. 3: The pressure-volume diagram for shock (blue curve) 
and detonation (red curve) in a gas mixture. The scaled en¬ 
ergy release is ( 7 — l)Q/ 7 -PoVb = 9, with the initial volume per 
unit mass Vo = 6.33 m 3 /kg and initial pressure Po = 1.33- 10 4 
Pa. The green straight line shows transition from the shock 
point s to the combustion products d in the CJ detonation. 

the jump from the point O (the initial state) to the point 
s at the shock, and then the reaction development behind 
the shock along the green straight line from the point s 
to the point d (detonation products). It may be shown 
that the C J detonation speed is determined by the energy 
release in the combustion processP 

The internal structure of the reaction zone in the det¬ 
onation front (transition from the point s to the point d 
in Fig. [ 3 ]) is specified by reaction kinetics of a particular 
fuel mixture. Here we illustrate the detonation struc¬ 
ture for the simplified model of acetylene-air combustion 
described by a single one-step Arrhenius reaction 22 

r)n 

— = - Kpaexp(-E a /T ), (7) 

where K = 10 9 m 3 /kgs is a pre-exponential factor, 
E a /T 0 = 29.3 is the scaled activation energy. Then the 
stationary detonation structure is determined by Eq. ([ 7 ]) 
rewritten in the reference frame of the front as 

3d 

u—= Kpaexp(-E a /T). (8) 

By solving Eq. ([s]) together with Eqs. 0-@ we find 
the detonation front structure as shown in Fig. [4j For 
obtaining the plots we have used the scaled energy re¬ 
lease ( 7 — 1)Q/(7-PoVo) = 9, with the initial gas pressure 
Po = 0.627 MPa and temperature T 0 = 293 K. Then by 
using the conservation laws Eqs. 0-0 one can find 
the maximal pressure in the detonation front achieved 
at the shock as P s = 3.7 MPa and maximal temperature 
achieved in the detonation products T d = 2988 K; these 
two values have been used for scaling in Fig. [4] As we can 
see, Fig. [4] presents the detonation structure similar to 
the qualitative sketch of Fig. [2j although the induction 
length is exhibited not so strongly by the acetylene-air 
detonations. 
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z(|im) 


P + pu 2 — Pq + PoD 2 , (10) 

/ _ P u 2 \ dT 

pu £ + Qa H-h — — k —— = 

\ p 2 J dz 

PoD f £q + Q 4-f \ > (11) 

V Po 2 J 

where e is the thermal energy per molecule, Q stands 
for Zeeman energy and a is the fraction of molecules in 
the metastable state. Zeeman energy Q, together with 
the activation energy E a , are determined by the applied 
magnetic field, and can be obtained by using Hamiltonian 
for Mni 2 -acetate molecule— 


FIG. 4: Characteristic detonation structure in a gaseous 
mixture with parameters representing acetylene-air combus¬ 
tion. The scaling is performed by using pressure at the shock 
P s = 3.7 MPa and the temperature in the detonation products 
Td = 2988 K; other computational parameters are specified in 
the text. 


III. MAGNETIC DETONATION STRUCTURE 
DUE TO HEAT TRANSFER 


A. Basic equations 


The theory of combustion detonation summarized 
briefly in the previous section provides the clue for de¬ 
scribing the magnetic detonation in crystals of nanomag¬ 
nets; conceptually the same model has been employed 
for the analysis of magnetic detonation in Ref. [2]. Still, 
here we stress important difference in some combustion 
and magnetic detonation properties, which is taken into 
account below. First, magnetic detonation in crystals of 
nanomagnets is extremely weak, so that it propagates 
with speed only slightly exceeding the sound speed cq in 
the crystal, D « cq. For that reason, in stark contrast 
to combustion detonations, the leading shock in magnetic 
detonation is not a discontinuity, but a smooth transition 
region with structure determined by transport processes, 
namely, thermal conduction k and volume viscosity 77 , 
which requires proper modifications of the conservation 
laws Eqs. 0-0 at the detonation front. 

So far the theoretical models of spin-avalanches in crys¬ 
tals of nanomagnets have taken into account only ther¬ 
mal conduction and neglected volume viscosityf^HlII in 
the present section we use the same approach and con¬ 
sider the influence of thermal conduction only. Volume 
viscosity will be taken into account in the next section. 
Although we deal with solid state, propagation of shocks 
and detonations in crystals of nanomagnets is also de¬ 
scribed by the hydrodynamic conservation laws of mass, 
momentum and energy similar to Eqs. 0-§ , see Refs. 

mm 


(9) 


H = -DS 2 -gg B B z S z , (12) 


where B z is the magnetic field, g = 1.94 is the gyromag- 
netic factor, p B is the Bohr magneton, and D = 0.65K 
is the magnetic anisotropy constant. We assume that 
the field is applied along the easy axis (z-axis). Zeeman 
energy is found as the energy difference between the sta¬ 
ble, S z = 10, and metastable, S z = —10, states of the 
nanomagnet as illustrated in Fig. |T] 

Q = 2gp B B z S^. (13) 

Zeeman energy increases linearly with the magnetic field. 
Activation energy is determined by the energy barrier 
between stable and metastable states, see Fig. [T] and for 
nanomagnets it may be found as 

E a = DS 2 - gp B B z S + ^pIB 2 z - (14) 

the activation energy is provided in temperature units. 
Dependence of Zeeman and activation energies on the 
applied magnetic field is presented in Fig. [5] for Mni 2 - 
acetate. As we can see, the activation energy decreases 
with the magnetic field and turns to zero at B z ss 10T; 
at higher fields the potential barrier disappears and the 
transition from the state S z = —10 to the state S z = 10 
is not hindered any more. 

It is convenient to define new dimensionless value 
r = p/po = Vo/V, which represents possible compres¬ 
sion of the crystal in the process of magnetic detonation. 
Weak shock and detonation cause only elastic deforma¬ 
tions, and the equation of state for the crystals of nano¬ 
magnets may be writt en as a combination of elastic and 
thermal components aJ^EH 


P 

Po 



n 


RATk B T a+1 r 
’ + M {a + 1) 0“ ’ 


(15) 


where M is the molecular mass with M = 1868 g/mol 
for Mni 2 -acetate, see Ref. [J. Similar to Ref. [2] we 
take n = 4, the Gruneisen coefficient T = 2, the problem 
dimension a = 3, and the Debye temperature 0 B = 38 
K. The coefficient A = 127r 4 /5 corresponds to the Debye 


pu = poD , 
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FIG. 5: The Zeeman energy Q and activation E a , both in 
temperature units, vs the external magnetic field B z . 


crystal lattice, and Co ~ 2000 m/s is sound speed in the 
crystal. Thermod ynam ic energy per one molecule of the 
crystal is given by^EH 


£= c|/r 
n 


-1 1 \ RAT a+1 

n- 1 + r ~ 7 + M(a + 1 ) 0 “ ' ^ ’ 


Equations (|9j)-( 11) provide a complete system for describ¬ 


ing magnetic detonation parameters in crystals of nano¬ 
magnets with thermal conduction taken into account. 

Solution to Eqs. Qd 11) is presented in Fig.[6]in terms 
of the pressure-volume diagram for the magnetic detona¬ 
tion structure similar to Fig. [3j We can see immediately 
that, since magnetic detonation is extremely weak, the 
shock curve (blue) and the detonation curve (red) almost 
coincide; similar to Fig. [3j the tangent line from the ini¬ 
tial state to the CJ detonation regime is shown by green. 
In the case of zero viscosity, pressure at the tangent line 
(label ”f”) follows from Eqs. ([9]) , (10) as 


Pt = P 0 + PoD 2 (1 - l/r). 


(17) 


In order to make difference between the shock and det¬ 
onation visible, we modify the pressure-volume diagram 
by extracting pressure at the tangent line, Eq. i.e. 

plotting P — P t versus V, see Fig. [7] Within such a 
representation, the tangent line (by definition) becomes 
simply the zero-line, while the shock and the detonation 
states may be well-distinguished and represented by two 
curves resembling parabola pieces. 


B. The analytical theory for the key magnetic 
detonation parameters 


As we can see in Fig. [TJ the limit of weak compression 
(r-l)Cl holds with a very good accuracy for magnetic 
detonations, and hence allows developing the analytical 
theory of the process. The system Eqs. ~([IT|) may be 
reduced to a single equation, which is equivalent to the 



FIG. 6: The pressure-volume diagram for the shock (blue) 
and magnetic detonation (red) in a Mni 2 -acetate crystal. The 
green line represents the tangent line from the initial state to 
the CJ regime; the magnetic field is B z = 4T. 



FIG. 7: The modified pressure-volume diagram for P — Pt for 
the shock (blue) and magnetic detonation (red) in a Mni 2 - 
acetate crystal. The reference pressure Pt of the green line 
corresponds to the tangent line Eq. (171 from the initial state 
to the CJ regime; the magnetic field is B z = 4T. 


Hugoniot relation 21 -^ with thermal conduction taken into 
account 


£ — Eq — 


1 r - 1 P 0 + P 

2 r p 0 


+ Q(& 


1) 


k dT 
PqD dz 


(18) 


By combining the equations of state for the crystals of 
nanomagnets, Eqs. (15), (16), we get rid of the ther¬ 
mal terms and obtain the expression for thermodynamic 
energy as function of P and r only 


1 P 
rT po 



nr 





n 

n — 1 


(r- 1) 


(19) 
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Then substituting Eq. (21) into the Hugonoit relation 


Eq. (18) we obtain pressure as 

p _ 2 r 

Po 2 — T(1 — r) 


{rQ(l - a) 


r — 1, 


£o 


k dT 
PqD dz 


n — 1 


(r ~ 1) - (r n - 1) 


n — 1 


( 20 ) 


In a similar way, by substituting Eqs. (15), (16) into Eq. 
(18) we find temperature 


rj- iQ; + l _ 


2{a + 1 )MQ% 
RA [2 - T(1 - r) 


{Q(l - a) 


n — 1 


(r - 1) - 


1 


r — 1 


r — 1. 


2 r 


-T + 1 Cq + 


n — 1 
k dT 
PqD dz 


( 21 ) 


Now we take into account that density variations in 
the magnetic detonation are small by introducing a small 
value S = r — 1 = p/po — 1 « 1. We can also neglect 
initial thermodynamic energy since Eq/Q ~ 10~ 6 . First, 
we analyze the detonation products, 2 —> oo, for which 
the transformation process is completed, dT/dz —> 0. 
Expanding Eq. (20) in powers of 6 and taking a = 0 


we obtain pressure in the detonation products as 

S d 


— =fqj 

Po 

1 

+ 2 




1 


^rQ ( 7T + 1 ) + Coin ~ 1) 


Si 


( 22 ) 


The scaled density deviation S d corresponding to the CJ 
point of the detonation products on Fig. [7]is unknown so 
far and we have to find it. At the CJ point the tangent 
line touches the detonation curve, which implies 


Pt =Pd , 


(23) 


(24) 


dPt^ _ dP d 

dr dr 

For small compression case <5 <C 1 we rewrite Eq. as 

Pt ~ Po 


Po 


= D 2 6{l-5), 


(25) 


Substituting Eqs.(22) and (25) intoEq.(23), (24), we find 


D 2 5 d (l-6 d ) =TQ + 


cq + rQ — + i 


i 

+ 2 


S d 


~ 2 ~Q f ^ + 1 ) + c o( n — 1) 


Si 


(26) 


D 2 (l-2S d ) =c 2 + TQ[- + l 


Q ) + c o( n — i) 


S d . (27) 


Then, eliminating D 2 from Eqs. (|26|, (|27|, and keeping 
terms as small as oc S d we obtain 

TQ (1 - 2TS d ) = 


+S 2 d 


rQ i + 


-5 


1 


c oi n + 1) 
(28)' 


We may also use the condition 4T 2 Q 5 d (n + l)cg jus¬ 
tified below to simplify Eq. (28) and find the density 


deviations in the detonation products as 


S d = - 
c 0 


1 / 2 rQ 


71+1 


(29) 


As an example, for the external magnetic field of B z = 4T 
we obtain the density deviations from the initial value 
as small as S d ~ 9.6 • 10 -3 <C 1. Temperature of the 
detonation products may be calculated with the accuracy 
of the first order terms in QS d as 


rjia-\-l _ 

j- A - 


(a + 1 )MQ% 


*9(1 + ijS d I + c 0 


2 n + 1 r 3 

o d 


12 


RA 

(a + 1 )QM0“ 
RA 


Figure [8] shows density, temperature and pressure of the 
detonation products versus the applied magnetic field as 
found in the numerical solution and predicted by the ana¬ 
lytical formulas Eq. (22), (291, (301; we observe excellent 



agreement of the theory with the numerical solution. We 
also point out that the magnetic detonation speed D 2 
may be obtained from Eq. (27) up to the leading terms 
in S d as: 


D 2 = 


1 


l-2S d 

T 2 


+ 


Cq + rQ ( — +1 

r 

2 


~ftQ ( tt + 1 ) + Coin ~ 1) 


(31) 


Keeping the leading terms we can rewrite Eq. (31) as: 
TQ 


D ~ Cq + 


2c q 


l + 1 ) + S i 


r Q/ r + i 
2c 2 U 


+ (n + 1) 
(32) 


Substituting Eq. (291 in Eq. (32) and taking into account 
TQ/c q Clwe find the magnetic detonation speed 


D 


co 


1 


TQ 


(33) 
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B(T) 



B(T) 


The tangent line intersects the shock curve at the point 
s, which implies P s = Pt or 


D 5 S (1 S s ) — 6 s c 0 


1+2^( n_1 ) 


(35) 


Taking into account the result for the magnetic detona¬ 
tion speed, Eq. (31), we obtain the equation for i5 s 


(1 — S s ) < Cq + 


+TQ - + 1 = (1 - 2 6 d ) 


~ttQ ( 77 + 1 ) + c 0 (n - 1) 


1 


S d 


i+ 7 $s{n- 1) 


, (36) 


which may be solved with the accuracy of TQ <C Cgi5 as 
5 S « 28 d (37) 


It is noted that density deviation just behind the discon¬ 
tinuous shock is about twice larger than in the detonation 
products, which may be seen already in Fig. [7 Then we 
find temperature on the shock similar to Eq. (30) 


+i _ (a + l)M0“ n + l^ §3 


rj-tCt- 1-1 _ 


RA 


12 


4(q + 1)TQMQ% I 2TQ 
3 cqRA 


(38) 


Density, temperature and pressure at the leading shock 
in magnetic detonation are shown in Fig. [8j we again 
observe a very good agreement of the analytical theory 
and the numerical solution. 


C. Isothermal discontinuity in magnetic detonation 


FIG. 8: Density, temperature and pressure at the leading 
shock and in the magnetic detonation products versus the ex¬ 
ternal magnetic field. Solid lines correspond to the numerical 
solution and the dashed lines represent the analytical theory. 


Thus, in agreement with Ref. [2], the magnetic detona¬ 
tion speed only slightly exceeds the sound speed in the 
crystals. 

It is also useful finding thermodynamic parameters of 
the crystal just behind the shock, a = 1, neglecting ther¬ 
mal conduction, i.e. treating the shock as a discontinuity 
similar to Ref. [2]. We will show in the next subsec¬ 
tion, that the same values for temperature, pressure and 
density hold at the isothermal discontinuity taking into 
account thermal conduction. So, we are looking for the 
density deviations 8 S corresponding to the crystal just 
behind an infinitely thin shock. Equation (20) yields up 
to the second order terms in <5 


Ps _ e 2 
— Sc 0 
Po 


1+ \ 6 (n~ !) 


(34) 


In this subsection we obtain the magnetic detonation 
structure in presence of thermal conduction. By using 
Eqs. @ we find temperature as a function of 

scaled density r inside the magnetic detonation front 


r+1= (« + P0‘ D 


TRAr 


D \ 1 - 


1 


- 7° ( r n - 1) 
n 


(39) 

By plotting temperature T according to Eq. ( f30j ) from 
the initial point O to the shock point s, we observe non¬ 
monotonic temperature variations, see Fig. [9j with the 
maximum value attained at the detonation point d. In 
that case, assuming hypothetic temperature variations 
along the red curve (solid and dashed) in Fig. [9] from 
O to d in the whole detonation front, one moves first 
from O to s ', then passes d to the shock point s, and 
then returns back from s to the final state d. Along such 
a pass we obtain the region with temperature decrease 
dT/dx < 0 on the way from d to the shock point s. 
However, the left hand-side of Eq. (18) turns to 0 only 
at the boundaries of the detonation wave z —> ±oo, i.e. 
it has to be either positive or negative in the transition 
zone from 0 to d , see Refs. I2II2B1 . Since detonation 
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FIG. 9: The isothermal discontinuity diagram. Temperature 
dependence on scaled density r according to Eq.( 391. Starting 
from the initial point O, all parameters vary continuously to 
s'. Then pressure and density experience isothermal disconti¬ 
nuity from the point s' to the point s (the black dashed line). 
The red dashed line corresponds to the nonphysical region in 
the shock phase. From the point s, the system relaxes to the 
final state d in the process of spin-flipping (magnetic detona¬ 
tion). Crystal compression at the point s' is r = 1.00012. 


definitely leads to temperature increase as compared to 
the initial state, then we have the condition dT/dx > 0 
satisfied everywhere in the magnetic detonation front, so 
that the hypothetic region with dT/dx < 0 is not physical 
and should be avoided. Then, variations of all values in 
magnetic detonation in Fig. [9] have to correspond to con¬ 
tinuous movement from the initial point O to the point 
s', which is followed by jump to the shock point s with 
subsequent continuous evolution to the final detonation 
point d. At the same time, because of thermal conduction 
in Eqs. temperature has to be continuous in 

magnetic detonation including the shock wave. For that 
reason, instead of a traditional shock, we obtain isother¬ 
mal discontinuity in magnetic detonation in the case of 
negligible volume viscosity when the shock structure is 
supported by thermal conduction only. An important 
point is that compression of the crystal in the point s' is 
quite small, r = 1.00012, so that the isothermal discon¬ 
tinuity resembles a shock for density and pressure. The 
visual difference is found only for temperature variations. 

The internal structure of magnetic detonation can be 
obtained by integrating the equation for kinetics of spin- 
relaxationp which resembles strongly the Arrhenius law 
of chemical kinetics Eq. Q, 

da a , „ , . , 

-Qj. = “-exp (~E a /T), (40) 


By adopting the reference frame of the stationary mag¬ 
netic detonation front, Eq. (40) can be written as 


= - -exp (~E a /T). 
ox T 


(41) 


By using Eq. (39), we solve Eq.(41) numerically starting 


FIG. 10: Magnetic detonation structure in Mni 2 -acetate 
crystals in the external magnetic field H = 4T. The plots 
are scaled by the temperature of the detonation products 
Td = 17.25 K and the shock pressure P s « 80kPa. Other 
parameters are specified in the text. 


from the initial point O to the point of isothermal dis¬ 
continuity s'. After that new values for density and pres¬ 
sure at the point s are calculated according to the jump 
conditions, with temperature and fraction of metastable 
molecules a kept constant. Then numerical calculation is 
performed from the point s to the final detonation state 
d. The result of the numerical solution is presented in 
Fig. [lO] for the Mni 2 -acetate crystals with the character¬ 
istic thickness of the detonation front L 0 = Cqt ~ 0.2 mm 
employed as the length scale. In Fig. 10 for illustra¬ 
tive purposes we plot the numerical solution for the 
thermal conduction coefficient k = 10 4 m 2 /s, which is 
about five times larger than the commonly accepted value 
k ~ 2-10 -5 m 2 /s, e.g. see Refs. [15I18| . Still, even in that 
case thickness of the heating region is much smaller than 
the region of spin-flipping, and fraction of the molecules 
in the metastable state is close to unity at the isothermal 
discontinuity, a = 0.993. 

Thus, as an intermediate conclusion, thermal conduc¬ 
tion influences only the temperature profile in magnetic 
detonation, with a minor effect on density and pressure, 
and with negligible modifications of the total front thick¬ 
ness. Much more dramatic modifications of the magnetic 
detonation structure are expected because of volume vis¬ 
cosity as shown in the next section. 


IV. DETONATION FRONT WITH VISCOSITY 

By taking into account volume viscosity, 77, we re-write 
the equations of mass, momentum and energy conserva¬ 
tion for magnetic detonation as 


Pouo = pu, 


(42) 


_ 2 2 du 

Po + p 0 u 0 =P + pu -77—, 


(43) 
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Po u o 



Po 

Po 



f p 

pu\ e H- 

V P 


u 2 \ dT du 

— — K— - VU-. 

2 / dx dx 


(44) 


We point out that, unlike commonly known shear viscos¬ 
ity arising due to relative motion of gas or fluid layers, 
volume viscosity describes momentum and energy dissi¬ 
pations due to compression of a medium. While shear 
viscosity is not typical for solid state processes such as 
magnetic detonation in crystals of nanomagnets, volume 
viscosity has to be considered. It is convenient character¬ 
izing the role of viscosity by the dimensionless parameter 


'/ - T 2 ' \ u / 

PqCqLq Po c o t 

This dimensionless viscosity plays a role conceptually 
similar to that of the inverse Reynolds number in fluid 
mechanicsP3 In particular, the parameter values rf = 
0.0025; 0.005; 0.01 employed below correspond to the do¬ 
main of the Reynolds numbers Re = 100 — 400. At such 
values of the Reynolds number gas and fluid flows are 
typically laminar; negligible role of viscosity is qualita¬ 
tively indicated by transition to turbulence, which hap¬ 
pens usually at larger values, Re ~ 10 3 and above. To the 
best of our knowledge, there have been no works, either 
experimental or theoretical, investigating volume viscos¬ 
ity in crystals of nanomagnets, and therefore we will take 
rf as a free parameter. 

We stress that even small values of volume viscosity 
rf ~ 0.01 lead to considerable changes in the magnetic 
detonation structure as shown in Figs. [TT] - [l4j Even on 
the pressure-volume diagram, Fig. (plotted, as before, 
for P — P t by extracting the CJ-tangent line for illustra¬ 
tive purposes) we observe that all discontinuous jumps of 
the previous modeP are replaced by continuous transition 
lines from the initial point to the CJ detonation prod¬ 
ucts. For very small viscosity, rf = 0.0025, the transition 
line (blue), although continuous, goes pretty close to the 
Hugoniot curve for the shock and then to the tangent line 
for the spin-flipping process. In a similar way, in Fig. |T3] 
the plot for rf = 0.0025 demonstrates quite abrupt ini¬ 
tial pressure increase - almost jump, which resembles 
strongly a discontinuous shock wave; this pressure in¬ 
crease is followed by pressure relaxation to a smaller 
value in the process of spin-flipping with Zeeman energy 
release. As we take larger values of volume viscosity, 
rf = 0.005; 0.01, deviation of the pressure-volume plots 
from the discontinuous shock model becomes much more 
pronounced, Fig. El and we find pressure and density 
in the shock wave decreasing as compared to the discon¬ 
tinuous case, Fig. [13] In fact, the very definition of a 
shock as a part of magnetic detonation front becomes 
ambiguous when volume viscosity is taken into account. 
To avoid the ambiguity, we notice that within the model 
of a discontinuous shock, pressure maximum is attained 
at the shock front, and then pressure goes down as the 
spin-flipping starts. In the same way it seems natural 



v/v 0 

FIG. 11: The pressure-volume diagram for different values 
of the scaled viscosity rf = 0.0025; 0.005; 0.01 in the external 
magnetic field B z = 4T. 



0 10 20 30 40 

z/L 0 


FIG. 12: Profiles of scaled pressure, temperature and fraction 
of nanomagnets in the metastable state for magnetic detona¬ 
tion for scaled viscosity rf = 0.005 in the external magnetic 
field B z = 4T. 


treating the point of maximal pressure in Figs. 12 13 as 
the back-side of the shock region. In Fig. [l3]we placed 
all three plots by choosing 2 = 0 as the position of the 
pressure maximum and hence the back side of the shocks. 
Then the regions corresponding to 2 < 0 belong to the 
shocks smoothed by viscosity, while the domain of z > 0 
may be treated roughly as the regions of spin-flipping 
with Zeeman energy release. Of course, such sep¬ 
aration is rather qualitative than quantitative, because 
in presence of considerable volume viscosity spin-flipping 
starts already in the shock, see Fig. [I2| corresponding to 
rf = 0.005, B z = 4T. In Fig. 12 the fraction of nanomag¬ 
nets in the metastable state is about a ss 0.6 at the point 
of pressure maximum. Figure |TT| compares temperature 
profiles for the scaled viscosity rf = 0.0025; 0.005; 0.01 
with the respective pressure maxima still placed at z = 0. 
As we can see, larger volume viscosity increases strongly 
temperature at the shock; in the case of rf = 0.01 tem¬ 
perature at z = 0 is only 5% smaller than the final det- 
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z/L 0 


FIG. 13: Pressure profiles in magnetic detonation for differ¬ 
ent values of the scaled viscosity rj' = 0.0025; 0.005; 0.01 in 
the external magnetic field B z = 4T. Pressure maxima are 
placed at the position 2 = 0, which may be defined as the 
back side of the shock region in the detonation. 


to observe steady well-developed magnetic detonation in 
the commonly experimental conditions. Instead, we sug¬ 
gest that most of the experimental points reported in Ref. 
PQ for ultra-fast magnetic avalanches correspond to mag¬ 
netic detonation in the process of development, which 
is also indicated by the average avalanche speed in the 
samples noticeably below the sound speed. Another im¬ 
portant feature of the experimental observations of Ref. 
ffl is that the ultra-fast avalanches were obtained for the 
magnetic field close to the quantum resonance values of 
the nanomagnets. Quantum resonances lead to strong 
decrease of the factor r in the kinetic equation of spin¬ 
flipping, Eq. (40), and hence of the characteristic width 
of the magnetic detonation front oc Lq = cqt, which may 
allow experimental observation of magnetic detonation 
in samples of conventionally employed sizes. 


V. SUMMARY 



z/L 0 


FIG. 14: Temperature profiles in magnetic detonation for 
different values of the scaled viscosity rf = 0.0025; 0.005; 0.01 
in the external magnetic field B z — 4T. Similar to Fig. |13| 
position 2 = 0 corresponds to the maximum pressure for each 
profile treated as the back side of the shock wave. 

onation temperature, i.e. it is about ~ 0.95T^. Tem¬ 
perature increase reduces strongly the size of the ’’pure” 
spin-flipping region, still, as we can see from Fig. [14J it 
is accompanied by strong increase of the width of the 
shock wave, so that total width of magnetic detonation 
front does not change much, being about ~ 10Lo for all 
three cases, i.e. about 1 — 2 mm in dimensional units. 

We point out that the experimentally employed sam¬ 
ple sizes for the crystals of nanomagnets are also about 2 
mm, e.g. see Refs. [Id- As a result, it is rather difficult 


In the present paper we have investigated the inter¬ 
nal structure of magnetic detonation in crystals of nano¬ 
magnets. Magnetic detonation is weak and propagates 
with speed only slightly exceeding the sound speed in 
the crystals. For that reason, in stark contrast to usual 
combustion detonations, transport processes - thermal 
conduction and volume viscosity - play an important role 
in forming the magnetic detonation structure. We show 
that, in the case of negligible volume viscosity, thermal 
conduction produces isothermal discontinuity instead of 
the leading shock in magnetic detonation. In the isother¬ 
mal discontinuity temperature of the crystal is continu¬ 
ous, while density and pressure experience jump. 

Volume viscosity leads to much more dramatic changes 
of the magnetic detonation structure as compared to the 
model of Ref. [2] with neglected transport processes. 
In that case all important thermodynamic parameters 
of the crystals acquire smooth profiles all over the mag¬ 
netic detonation front including the leading shock. In 
addition, the very concept of the leading shock requires 
unambiguous definition - here we suggest to specify the 
back side of the shock as the position of the pressure max¬ 
imum. As the relative role of volume viscosity increases, 
the leading shock becomes wider and may exceed consid¬ 
erably the zone of spin-flipping by size. Still, total size of 
a magnetic detonation front does not change much with 
variations of volume viscosity since decrease of the spin¬ 
flipping region is compensated by increase of the shock 
width. 
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